Evidence for 2kp Electron-Electron Scattering Processes in Coulomb Drag 
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Measurements and calculations of Coulomb drag between two low density, closely spaced, two- 
dimensional electron systems are reported. The experimentally measured drag exceeds that calcu- 
lated in the random phase approximation by a significant, and density dependent, factor. Studies 
of the dependence of the measured drag on the difference in density between the two layers clearly 
demonstrate that previously ignored q — 2kF scattering processes can be very important to the drag 
at low densities and small layer separations. 
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Interlayer Coulomb interactions between electrons in 
spatially separated two-dimensional electron systems 
(2DES) can lead to the condensation of exotic bilayer 
collective phases which do not occur in single layer sys- 
tems [0. A dramatic example of this occurs when the 
total number of electrons in the bilayer system equals 
the number of states in the lowest spin-resolved Landau 
level created by a perpendicular magnetic field. If the 
separation between the layers is small enough, intra- and 
interlayer Coulomb interactions can combine to produce 
a spontaneous interlayer phase coherent state in which 
there is complete quantum uncertainty as to which layer 
any given electron is in. This unusual state has been 
predicted to possess a number of remarkable proper- 
ties including counter-flow superfluidity and Josephson- 
like interlayer tunneling. Recent experiments ||^,^ have 
strongly supported these predictions. 

In spite of the key role played by interlayer electron- 
electron interactions in the stabilization of bilayer con- 
densed phases, there have been very few quantitative 
measurements of their strength. One relatively new tech- 
nique [^j5| for making such measurements consists of 
recording the "drag" voltage Vd which develops in one 
electron layer in response to a current flow I confined 
solely to the other layer. The resulting drag resistance 
Rd = Vd/I is directly proportional to the interlayer 
momentum relaxation rate in the system. At low tem- 
peratures and with closely spaced layers, direct inter- 
layer Coulomb scattering dominates this rate [P^, but 
in general other processes, such as virtual phonon ex- 
change ||7|-|Tot and plasmon-enhanced Coulomb scattering 
pl|-|r3t, can also contribute. 

At temperatures T small compared to the Fermi tem- 
perature Tp the divergent phase space for forward {q = 0) 
and backward [q = 2kF) Coulomb scattering in a clean 
2DES leads to ln{T) corrections to the usual depen- 
dence of the inverse thermal quasiparticle lifetime [ p^ . 
In drag, however, the situation is somewhat different. 
Most importantly, the Fourier transformed bare inter- 
layer Coulomb interaction is exponentially sensitive to 
the spacing d between the layers: V{q) ^ e^'^'^/q. This 



effectively suppresses scattering processes with momen- 
tum transfers in excess of q ^ 1/d. Thus, if kpd >> 1 
backward scattering q = 2kp processes are unimpor- 
tant. In this case the drag is dominated by small an- 
gle events and one expects Rd ^T"^- This situation 
was appropriate in most early Coulomb drag experiments 
P,P|JT^. In this case the IniT) dependence is absent be- 
cause the momentum relaxation rate is determined by 
the electron-electron scattering cross-section weighted by 
1 — cos9 q"^, where 9 is the scattering angle. This factor 
vanishes rapidly enough at small q to suppress the low 
angle phase space divergence. 

Backward scattering processes can become important 
in drag at low carrier densities and for closely spaced 
layers since then kpd may be of order or even less than 
unity. In this case a T^ln(T) temperature dependence 
of the drag resistance is expected at low T. Although 
definitive identification of such logarithmic corrections 
is difficult, we nonetheless report here strong evidence 
for q = 2kp Coulomb drag scattering processes in low 
density bilayer systems. This evidence is obtained not 
from the temperature dependence of the drag, but rather 
from its sensitivity to the different densities in the two 
2D layers. 

The drag experiments reported here were performed on 
bilayer 2DESs in GaAs/AlGaAs heterostructures. These 
samples contain two 18nm GaAs quantum wells sepa- 
rated by a lOnm Alo.gGao.iAs barrier layer. This double 
well structure is sandwiched between thick Alo.3Gao.7As 
layers which contain Si (5-doping layers positioned about 
220nm from the GaAs quantum wells. At low tempera- 
tures each quantum well contains a 2DES of nominal den- 
sity 5.2 X 10^°cm~^ and mobihty 1 x lO^cm^/Vs. Stan- 
dard photolithographic techniques were used to pattern a 
mesa, 40/im wide by 400/im long, on the sample. Diffused 
In ohmic contacts were placed at the ends of arms which 
extend outward from this bar-shaped central mesa. A se- 
lective depletion scheme allows each of these contacts to 
be connected to the central region through either 2D layer 
separately Metal gate electrodes on both sides of 
the thinned hctcrostructure sample provide control over 
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FIG. 1. Drag resistivity vs. temperature for six different 
densities A'o- Densities are in units of lO^^cm"^. Dashed 
line is a least-squares fit of the A^o = 3.1 x lO^^cm"^ data to 
Pd = AT"". 



the electron densities A''i and N2 of each 2D layer. Drag 
measurements were performed by driving a current, typ- 
ically lOnA at 13Hz, down the bar through one 2D layer 
while the drag voltage which develops in the other layer 
is recorded. Considerable care was exercised in order to 
eliminate spurious contributions to the drag signal aris- 
ing from the finite tunneling resistance (> lOOMSl) and 
capacitance 140pF) between the layers. No effect on 
the drag resistance was found when the role of drive and 
drag layers were interchanged. 

Figure 1 shows the measured drag resistivity po vs. 
temperature at six different balanced (i.e. Ni — N2 = 
No) densities in the two 2D layers: Nq = 1.7, 2.3, 3.1, 3.8, 
4.7, and 5.2 x 10^°cm~^. As expected, the drag resistance 
increases with temperature and is larger at lower density. 
The dashed line in the figure represents an unweighted 
least-squares fit of the iVg = 3.1 x lO^^cm"^ drag data to 
a simple quadratic temperature dependence: pu — AT^. 
The fitted coefficient A exceeds the simplest theoretical 
estimate of Coulomb drag by a factor of 5.9 at this den- 
sity. This estimate is based upon a model which as- 
sumes two ideally thin 2D layers separated by the present 
sample's center-to-center spacing oi d = 28nm, low tem- 
peratures {T << Tp), and a predominance of small- 
angle scattering {kpd >> 1). The model also treats 
screening of the interlayer Coulomb interaction in the 
random phase approximation (RPA) under the assump- 
tion QTpd >> 1, with qxF the Thomas-Fermi screening 
wavevector. We find that this model underestimates the 
drag at all densities studied and that the shortfall in- 
creases from about a factor of 2 at A^o = 8.8 x lO^^cm^^ 
to a factor of 10 at A'^o = 1-7 x lO^^cm^^^ 

The discrepancy between theory and experiment in the 
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FIG. 2. Drag resistivity vs. density at three temperatures: 
T = 4K, 2K, and IK. Solid lines are proportional to Nq'^; 



dashed line proportional to A'q 



magnitude of the drag is substantial but not unusual. 
Similar discrepancies have been reported in electron- 
electron Q, electron- hole and low density hole- hole 
p6[ samples. Although phonon exchange can contribute 
to the drag between 2D electron systems [7|-pX[|, it is un- 
likely to be very important here. Experiments in which 
the phonon contribution was deemed to be comparable to 
the Coulomb scattering component show a much smaller 
drag resistance than we report here [f 10|. Detailed cal- 
culations support this conclusion Furthermore, the 
phonon contribution has been traditionally identified via 
the peak it produces in p^/T'^ vs. T. This peak, which 
occurs when the mean thermal phonon wavevector is 
comparable to 2k p, is not observable in the data pre- 
sented here. 

Another possible source of the enhanced drag which we 
observe are higher order many-body effects not captured 
by an RPA treatment of screening. Such effects generally 
become more important at low density, but there have 
been few quantitative estimates of their importance in 
Coulomb drag |]l^. To investigate this further we show 
in Fig. 2 a log-log plot of the drag resistivity ppi vs. 
density A^o at three fixed temperatures: T = 1,2, and 4K. 
The straight solid lines, while merely guides to the eye, 
represent a Nq'^ dependence. In contrast, the straight 
dashed line shows the N^^ dependence expected from the 



simple model described above y,|6[. It is apparent that 
our experimental results are much better approximated 
by the quartic density dependence. 

It is also clear from Fig. 1 that the drag data are not 
well fit by the temperature dependence of the simplest 
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FIG. 3. Sensitivity of drag to antisymmetric density 
changes AA'' in the two 2D layers. Left panels: Average den- 
sity iVo = 3.7xl0^'^cm"^. Right panels: A^o = 5.2 x 10^°cm"^. 
For each density, data from two temperatures is shown. Inset: 
Density dependence of the cross-over temperature Tc- 

theoretical model. Indeed, we find that an excellent fit 
to the A^o = 3.1 X 10-^°cm~^ data is provided by pn — 
AT'^ln{T/To) with A = -OAUQ/a and Tq = 1160K. 
We emphasize, however, that other fitting functions also 
work fairly well {e.g. po = BT^'^) and that we are not 
proposing any specific analytic form. 

We now turn to the dependence of Coulomb drag on 
antisymmetric changes in the density of the two layers: 
Ni^ Nq + an/2 and N2 = No - AN/2. Such changes 
are readily imposed by applying a small dc bias voltage 
between the two 2D layers Figure 3 illustrates the 
effect of such density imbalances on the drag in our sam- 
ples. Data at two different average densities Nq and two 
temperatures are shown. As the figure makes clear, very 
different behavior is observed at low and high tempera- 
tures. At high T the drag is found to increase, roughly 
quadratically, with AA^. In contrast, low temperatures 
produce the opposite result: the drag falls, again roughly 
quadratically, with AA^. We find a smooth transition be- 
tween the two regimes and a well-defined temperature Tc 
at which the drag is roughly independent of AA^ for small 
AA^. The inset to the figure suggests an approximately 
linear dependence of the cross-over temperature on 
the average layer density A^o- In terms of Fermi temper- 
atures, we find that Tc/Tp « 0.12 roughly defines the 
cross-over temperature for this sample. 

We propose that the results shown in Fig. 3 strongly 
suggest that g = 2/cf electron-electron scattering pro- 
cesses are important to the Coulomb drag in the present 



sample. To see this, we begin by noting that in the 
simple theoretical model presented earlier for compar- 
ative purposes, drag increases with AA^. This is a re- 
sult of the fact that the drag resistivity in that model is 
proportional to {NiN2)^^^'^ . For small AA^, the model 
leads to Apo/pn = +3iAN/Na)'^ /8. But the validity of 
this model depends, in part, upon the assumption that 
kpd >> 1 and the resultant restriction to small angle 
scattering processes. For the present sample, however, 
kpd ranges from about 0.9 to 1.6 for the densities used. 
These values are sufficiently small that 2k p scatterings 
cannot be ignored, especially in view of the large phase 
space for such events. Furthermore, back-scattering pro- 
cesses offer a natural way to understand the imbalance 
dependence of the drag at low temperatures. For drag 
it is the product of the phase space availability in each 
layer which matters. If the two 2D systems have the same 
density, their Fermi surfaces have the same diameter and 
the phase divergences at q — 2kp reinforce one another. 
When a density imbalance AA^ is imposed, the Fermi sur- 
faces no longer overlap and the joint phase space product 
is diminished, dramatically so at low temperature, and 
this causes the drag to decrease. 

At higher temperatures the phase space singularities 
are washed out. Owing to the weighting of the scat- 
tering cross-section, the effect of this thermal smearing 
is particularly important for large angle q — 2kp pro- 
cesses. Indeed, the relative importance of backward scat- 
tering events declines at high temperatures and the mean 
scattering angle diminishes. For these reasons, it seems 
plausible that the drag might increase with density im- 
balance AA^ in a way similar to that which occurs un- 
der the assumptions of the small-angle scattering theory 
discussed above. In order to make this qualitative argu- 
ment for a change in sign of the AA^ dependence of drag 
more compelling, we have performed detailed numerical 
calculations of the drag using the theoretical framework 
originally employed by Gramila, et al. Q and Jauho and 
Smith H. Ahhough our calculations follow the same 
Boltzmann equation approach as in these earlier 
works, we have applied them here to the significantly 
different regime of sample parameters (layer separation 
and density) appropriate to the present experiments. 

Figure 4 shows typical results of our calculations. In 
the main panel the drag at a density of A^o = 3.7 x 
10^"cm~^ is shown. These results include the effect of 
the finite thickness of the 2D systems and incorpo- 
rate screening at the RPA level. No attempt was made to 
include phonon or plasmon-related contributions to the 
drag. As expected, the magnitude of the calculated drag 
falls short of the experimental results at the same density 
by a significant factor 5). Nonetheless, as the insets 
demonstrate, the calculations do reproduce the change in 
sign of the dependence of drag on antisymmetric density 
changes: A^i = A^o + AN/2 and A^2 = A^o - AA^/2. In 
common with the experiment, at low temperatures po 
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FIG. 4. Calculated Coulomb drag vs. temperatures for 
A^o = 3.7 X lQ^^cm~^. Insets show dependence of drag on 
antisymmetric density changes AA'^ at low and high temper- 
atures. 



FIG. 5. Drag intensity h{q), divided by T^, vs. momentum 
transfer g at T = 1.4K and 9.0K, for No = 3.7 x 10^° cm~^ 
Note the strong peak near q — 2kF at low temperature. Drag 
resistivity is the area under these curves. 



falls roughly quadratically with AA^ while at high tem- 
peratures it rises. The calculated cross-over temperature 
Tc is about 6.3K; this is about a factor of 4 higher than 
the experimental value at the same density. It is quite 
clear from the calculations that the drag at low temper- 
atures contains a strong component from Coulomb back- 
scattering processes sharply peaked around q = 2kp. At 
high temperatures this component is reduced and a broad 
distribution of smaller scattering angles dominates the 
drag. This is clearly demonstrated in Fig. 5 where the 
drag "intensity" h{q) is plotted v.s. momentum transfer q. 
(The net drag, pz?, is obtained by integrating h{q) over all 
q.) Thus, these calculations strongly support the qualita- 
tive argument given above and demonstrate that q = 2kp 
backward scattering processes can be very important at 
low temperatures in samples with low electron densities 
and small layer separations. These processes are directly 
detectable via the unusual dependence of the drag on 
antisymmetric density changes in the double layer 2D 
system. 

In conclusion, we have measured Coulomb drag in low 
density 2D electron systems with small layer separations. 
Our results show that the drag is substantially larger 
than theoretical results based on RPA screening of the 
interlayer Coulomb interaction. We find the level of dis- 
agreement to grow steadily worse as the density is re- 
duced. At low temperatures our experimental data and 
numerical calculations clearly demonstrate that previ- 
ously ignored large angle q — 2kp scattering processes 
can be quite important in Coulomb drag. 
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